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ABSTRACT 

Since the discovery of the first double neutron star (DNS) system in 1975 by Hulse and Taylor, there are 
currently 8 confirmed DNS in our galaxy. For every system, the masses of both neutron stars, the orbital semi- 
major axis and eccentricity are measured, and proper motion is known for half of the systems. Using the orbital 
parameters and kinematic information, if available, as constraints for all system, we investigate the immediate 
progenitor mass of the second-born neutron star and the magnitude of the supernova kick it received at birth, 
with the primary goal to understand the core collapse mechanism leading to neutron star formation. Compared 
to earlier studies, we use a novel method to address the uncertainty related to the unknown radial velocity of the 
observed systems. For PSR B1534+12 and PSR B1913+16, the kick magnitudes are 150 - 270 km/s and 190 - 
450 km/s (with 95% confidence) respectively, and the progenitor masses of the 2nd born neutron stars are 1.3 - 
3.4 Mq and 1.4 - 5.0 M Q (95%), respectively. These suggest that the 2nd born neutron star was formed by an 
iron core collapse supernova in both systems. For PSR J0737-3039, on the other hand, the kick magnitude is 
only 5 - 120 km/s (95%), and the progenitor mass of the 2nd born neutron star is 1.3 - 1.9 M Q (95%). Because 
of the relatively low progenitor mass and kick magnitude, the formation of the 2nd born neutron star in PSR 
J0737-3039 is potentially connected to an electron capture supernova of a massive O-Ne-Mg white dwarf. 
For the remaining 5 Galactic DNS, the kick magnitude ranges from several tens to several hundreds of km/s, 
and the progenitor mass of the 2nd formed neutron star can be as low as ~1.5 M Q , or as high as ^8 M Q . 
Therefore in these systems, it is not clear which type of supernova is more likely to form the 2nd neutron star. 
Subject headings: binaries: close — pulsars: individual(PSR B1534+12, PSR B1913+16, PSR J0737-3039, 
PSR J 1 5 1 8+4904, PSR J 1 756-225 1 , PSR J 1 8 1 1 - 1 736, PSR J 1 829+2456, PSR J 1 906+0746) 
— stars: evolution — stars: kinematics — stars: neutron — supernovae: general 



1. INTRODUCTION 

Over the years, computation modeling and hydrodynamic 
calculations of supernova explosions due to core-collapse of 
massive stars have improved our theoretical understanding of 
how compact objects form. At present, two types of super- 
novae are actively discussed in the literature: (i) After iron 
is synthesized within a massive star, the core completes its 
final nuclear burning state. It then becomes degenerate and 
keeps growing in mass due to silicon shell burning outside 
the core. Eventually, it crosses the Chandrasekhar mass limit 
and ensues collapse. When the core reaches nuclear density, 
it bounces back and explodes the star as a supernova. (Bethe 
1990; Mezzacappa 2005; Woosley & Bloom 2006; Kotake et 
al 2006; Janka et al. 2007) (ii) On the other hand, if the core is 
not massive enough to ignite neon, nuclear burning within the 
core stops. The core becomes degenerate due to cooling by 
neutrino emission. Because of carbon shell burning outside 
the core, the mass of the core approaches the Chandrasekhar 
limit. When the core density and temperature rise above a cer- 
tain threshold, electron capture by nuclei (such as 20 Ne, 20 F, 
2A Na, and 2A Mg) takes place, decreasing degenerate electron 
pressure support. As a result, the core contracts rapidly and 
eventually bounces back, driving a supernova explosion. This 
scenario is referred to as electron capture supernova (ECS) 
(Miyaji et al. 1980; Nomoto 1987). Although the general idea 
of ECS has been established for several decades, the mecha- 
nism of how electron capture occurs in these relatively less 
massive cores is not yet settled. 

Regardless of the supernova mechanism, it is expected that 

TsingWong2012@u. northwestern.edu, b-willems@northwestem.edu, vick; 



some asymmetries develop during the collapse, imparting a 
recoil kick to the nascent neutron star (see Harrison et al. 
1993; Frail et al. 1994; Lyne & Loimer 1994). The observed 
proper motions of 233 Galactic pulsars show that there is a 
class of neutron stars receiving high recoil velocities at birth 
(Hobbs et al. 2005). On the other hand, Pfahl et al. (2002) 
studied the orbital parameters of low eccentric Be/X-ray bi- 
naries, and pointed out that these neutron stars could only 
have received low kick velocities (< 50 km/s) at the time of 
formation. There is thus evidence for the existence of two 
different neutron star formation mechanisms (van den Heuvel 
2004, 2007), potentially distinguishable by the magnitude of 
the natal kick of the neutron star. 

Currently, there are 8 DNS observed in our galaxy. Their 
orbital parameters are summarized in table 1 . The magnitude 
of the natal kick and the progenitor mass of the second-born 
neutron star (NS2) in these systems, except PSR J 1906+0746, 
have been studied previously by Willems et. al. (2004, 2006), 
Thorsett et. al. (2005), Stairs et al. (2006), and Wang et 
al. (2006). In their study, Wang et al. (2006) studied the 
formation of all DNS except PSR J 1906+0746, and derived 
constraints on the natal kick magnitude and progenitor mass 
of NS2 by studying the orbital dynamics of asymmetric su- 
pernova explosions, but not their kinematic history. The other 
authors focused on PSR B 1534+ 12, PSR B 19 13+ 16, and PSR 
J0737-3039 and used both orbital dynamics of asymmetric su- 
pernova explosions and the measured proper motion of these 
systems to constrain the natal kick magnitude and progenitor 
mass of NS2. The authors traced the motion of the systems 
in the Galaxy backwards in time and assumed each crossing 
of the binary through the Galactic disk within the age of the 
lrMveStetaeBa a possible birth site of the DNS. However, the 
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backward calculation of the motion in the Galaxy requires the 
knowledge of the radial velocity which cannot be measured. 
Willems et al. (2004, 2006), Thorsett et al (2005) and Stairs 
et al. (2006) therefore considered possible radial velocities 
drawn from either a uniform or Gaussian radial velocity dis- 
tribution. 

In this paper, we avoid the assumption of a present-day ra- 
dial velocity distribution by carrying out Monte Carlo simula- 
tions of the motion of the observed DNS in the Galaxy that is 
forward in time instead of backward. For this purpose, we dis- 
tribute populations of DNS progenitors in the Galaxy accord- 
ing to a double exponential distribution function. Each system 
in the populations is assigned a velocity equal to the vector 
sum of the local Galactic rotational velocity and an isotropic 
kick velocity with a magnitude generated from the kick ve- 
locity distribution function, derived from the supernova or- 
bital dynamics constraints for each of the observed systems. 
Compared to Willems et al.(2004), Thorsett et al. (2005), and 
Wang et al. (2006), we also relax the constraint that NS2's 
progenitor needs to be more massive than 2.1 M Q , which is 
a convectional limit for a helium star that explodes in a core 
collapse supernova. Instead, we restrict the progenitor to be 
more massive than NS2 only. Furthermore we study all eight 
known DNS in our galaxy, using the up-to-date observational 
parameters listed in Table 1 . 

The methodology of the adopted analysis is outlined in 
more detail in §2, while results for the individual systems 
and comparison with earlier studies are discussed in § 3. In 
§ 4, we summarize our results and discuss their implications 
for the supernova forming the second-born neutron star in the 
observed DNS. 

2. METHODS OF CALCULATION 

Tauris & van den Heuvel (2004) reviewed the general sce- 
nario for forming DNS. It starts with a binary consisting of 
two massive ZAMS stars ( > 10M Q ). The primary star, which 
is initially the more massive one, will leave the main sequence 
first and explode as the first supernova, in which it forms the 
first-born neutron star (NS1) in this binary. There may be 
mass transfer from the primary to the secondary star before 
that supernova takes place. As the secondary star leaves the 
main sequence, wind mass transfer from the secondary to NS 1 
can occur, leading to the formation of a X-ray binary. While 
NS 1 is accreting mass from its companion, it is spun up or "re- 
cycled", and its spin axis is expected to align with the orbital 
angular momentum. Later on, the X-ray binary may end up 
in a common envelope phase, during which NS 1 is engulfed 
by the extended envelope of its companion. The orbit shrinks 
and circularizes rapidly due to the frictional forces acting on 
NS 1 as it moves through the envelope of the companion. The 
lost orbital energy is deposited into the gas envelope, causing 
it to become gravitationally unbound and get ejected, leaving 
behind the companion's bare helium core. The helium star 
eventually explodes in a supernova and becomes NS2, orbit- 
ing the recycled, milli-second pulsar NS 1 in a compact binary 
orbit. Mass transfer from the helium star to NS 1 can also oc- 
cur prior to the second supernova explosion. 

Here, our goal is to use the current DNS binary properties 
listed in Table 1 to constrain as best possible the magnitude of 
the natal kick and the mass of the progenitor of NS2 in the 8 
known DNS in our galaxy at the time of the second supernova. 
We also predict the unknown radial velocity of each system, 
as well as the total transverse velocity and NS 1 spin tilt angle 
for those systems where these properties have not been mea- 



sured yet. We achieve these goals by means of Monte Carlo 
simulations, incorporating orbital dynamics during supernova 
explosions and kinematic of binaries in the Galactic potential, 
as detailed in what follows. 

We start our analysis by deriving a probability distribution 
function (PDF) for the natal kick of NS2, which depends on 
the semi-major axis and orbital eccentricity immediately af- 
ter the second supernova (A post sN and e pos sN respectively). As 
the tidal interaction between two neutron stars in DNS is ex- 
tremely weak (Bildsten & Cutler 1992), the major mecha- 
nism that changes the orbital semi-major axis and eccentric- 
ity of DNS is the emission of gravitational waves. To account 
for gravitational radiation driven orbital evolution, we use the 
current observational parameters as initial conditions and in- 
tegrate equations (35) and (36) in Junker & Schafer (1992) 
backwards in time. Because of the uncertainty in the true age 
of the observed DNS (see Kiziltan et al. 2009 for details), we 
randomly draw the integration time f, from a uniform distri- 
bution of ages between and r c , where t c is the characteristic 
age of the DNS. If r c is greater than the age of our galaxy, 
which is about 10 Gyr , we use 10 Gyr as the upper limit 
for h instead. For PSR B1534+12, PSR B1913+16, and PSR 
J0737-3039, different upper limits are adopted for f„ as dis- 
cussed in more detail in the next section. 

Once Ap OSt sN and e post sN are known from the gravitational 
radiation integration, we randomly draw the natal kick veloc- 
ity (Vk) from a uniform distribution of values between and 
2500 km/s (i.e., we assume no prior kick distribution derived 
from pulsar samples analysis), and the kick direction (8,4>) 
from an isotropic distribution. Here 9 is the polar angle be- 
tween the natal kick velocity and the instantaneous orbital ve- 
locity of NS2's progenitor at the time of the second super- 
nova explosion, and <p is the corresponding azimuthal angle 
(see Figure 1 in Kalogera 2000 for a graphical representation). 
Using the conservation laws of orbital energy and angular mo- 
mentum expressed by equations (19) and (20) in Willems et 
al. (2005), we obtain the orbital semi-major axis (A pre sN) im- 
mediately before the second supernova and the mass of NS2's 
helium star (M2;). The pre-supernova orbit is assumed to be 
tidally circularized at the time of the second supernova explo- 
sion. On the other hand, the spin tilt angle of NS1, which 
is the angle between NS 1 axis and the post supernova orbital 
angular momentum axis, can be calculated by following equa- 
tion (3) in Kalogera (2000). Also, using equations (12) in the 
same paper, we find the polar angle 8' between the natal kick 
and the pre-supernova orbital angular momental axis. 

Immediately before the second supernova explosion, there 
may be mass transfer from NS2's progenitor to NS1. If the 
radius R21 of NS2's progenitor is greater than its Roche lobe 
radius Ri, which is approximated by equation (2) in Eggleton 
(1983), Roche lobe overflow (RLO) takes place. Since NS2's 
progenitor is a bare helium star before the second supernova, 
we use equation (A10) in Kalogera & Webbink (1998) to cal- 
culate Rxi- Note that the left hand side of that equation should 
be logRfj e .f instead of Rfie.f- In order to have dynamically 
stable RLO, so that NS2's progenitor does not coalesce with 
NS 1 to form an isolated black hole, the mass ratio M2, to M\ 
needs to be less than 3.5 (Ivanova et al. 2003). Even if NS2's 
progenitor does not overflow its Roche lobe, there may be 
mass transfer by stellar wind. In addition, M2, is restricted to 
be less than 8 M Q , which is a conventional limit that a helium 
star will become a neutron star instead of a black hole (see, 
e.g., Figure 1 in Belczynski et al. 2002; Table 16.4 in Tauris 
& van den Heuvel 2004). 
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TABLE 1 

Parameters of the 8 known DNS in our galaxy 



System 



Ma 



Mi 



Mi ■ 



M 2 > 



e C ur 3 



PSRB1534+12 rl 

PSRB1913+16 r2 

PSR J0737-3039 13 

PSRJ1518+4904 14 

PSRJ1756-2251' 5 

PSRJ1811-1736' 6 

PSRJ1829+2456 r7 

PSRJ1906+0746 r8 



15 37 09.96 
19 15 28.00 
07 37 51.25 
15 18 16.80 

17 56 46.63 

18 11 55.03 

18 29 34.60 

19 06 48.67 



11 55 55.55 
16 06 27.40 
-30 39 40.71 
49 04 34.25 
-22 51 59.40 
-17 36 37.70 
24 56 19.00 
07 46 28.60 



1.02 

8.3(1.4) 
1.15 
0.625 
2.5 
6.0 
1.2 
5.4 



1.34(1) 
-3.27(35) 
-3.82(62) 
-0.67(4) 
-0.7(2) 



-25.05(2) 250 
-1.04(42) 110 
2.13(23) 210 
-8.53(4) 20000 
443 
1830 
■ ■ ■ 12400 
••• 0.112 



1.3332(10) 
1.4408(3) 
1.337(5) 
72+ ' 51 

1.312(17) 
! 62 +0.22 

1 - IH -0.48 
1.365(18) 



1.3452(10) 
1.3873(3) 
1.250(5) 
2 00 +0 58 

1 9C0+0.018 

136^ 
1.248(18) 



3.28 0.274 25(155) ± 3.8 

2.80 0.617 18(162) ±6 

1.26 0.0878 <15 1 

24.7 0.249 

2.70 0.181 

40.7 0.828 

6.36 0.139 

1.75 0.0853 



a Right ascension (J2000.0) (hours, minutes, and seconds) 
b Declination(J2000.0) (degrees, arcminutes, and arcseconds) 
c Distance (kpc) 

d Proper motion in R.A. (mas yr -1 ) 
e Proper motion in dec. (mas yr -1 ) 
f Characteristic age (Myr) 
8 Mass of 1st bornNS(M Q ) 
h Mass of 2nd born NS (M Q ) 
1 Current semimajor axis (Rq) 
' Current orbital eccentricity 

k Misalignment angle between pulsar spin axis and post-supernova orbital angular momentum axis (degree). The number in parenthesis correspond to another 
possible angle 

1 Different data analysis model gives a different upper limit on this angle. See section 4.3 for details. 
rl Wolszczan 1991; Arzoumanian et al. 1999; Stairs et al. 2002; Konacki et al. 2003; Stairs et al. 2004 

r2 Hulse & Taylor 1975; Taylor et al. 1976, 1979; Taylor & Weisberg 1982, 1989; Damour & Taylor 1991; Arzoumanian et al. 1999; Wex et al. 2000. 

r3 Burgay et al. 2003; Lyne et al. 2004.; Lorimer et al. 2007; Ferdman 2008 et al. 2008; Breton 2009; Deller et al. 2009 

r4 Thorsett & Chakrabarty 1999; Janssen et al. 2008 

15 Faulkner et al. 2005; Ferdman 2008 

r6 Lyne et al. 2000; Stairs 2004; Corongiu et al. 2007 

r7 Champion et al. 2004, 2005; Stairs 2004; 

r8 Lorimer et al. 2006, Kasian 2008, Stairs 2008 



The conservation laws of orbital energy and angular mo- 
mentum allow for real solutions only for M 2( -, Vk, 0, <f>, and 
ApreSN values satisfying the orbital dynamics constraints given 
by equation (21)-(27) in Willems et al. (2005). For PSR 
B1534+12, PSR B1913+16, and PSR J0737-3039, these con- 
straints are supplemented with the measured NS1 spin tilt 
angle constraints. Applying these constraints together with 
the constraints imposing stability of any mass transfer tak- 
ing place right before the second supernova explosion yields 
a PDF of the natal kick velocity of NS2. A corresponding 
PDF for the kick velocity Vk. sys imparted to the binary's center 
of mass is obtained using equations (28)-(32) in Willems et 
al. (2005). In all past studies, kinematic analyses have either 
been skipped (Wang et al. 2006) or simplified to one dimen- 
sion (Piran & Shaviv 2004) or have suffered by ad hoc as- 
sumptions about the un-measurable radial velocities. Willems 
et al. (2006) used general DNS population synthesis mod- 
els to guide this assumption. In this sutdy we employ a new 
method that still uses forward in time kinematic simulations, 
but they are tied to each observed system instead of being 
general. 

We perform a Monte Carlo simulation of the motion of 
DNS in the Galactic potential. For this purpose, we adopt 
a reference frame with origin at the Galactic center, with z- 
axis pointing to the northern Galactic pole, and with the x- 
axis pointing in the direction from the Sun to the Galactic 
center. With respect to this reference frame, the Sun is lo- 
cated aXR Q = (-8.5,0,0.03) kpc (Ghez et al. 2008, Gillessen 
et al. 2009, Reid et al. 2009, Joshi 2007) and has a pe- 
culiar velocity of (10,5,7)km/s (Bienayme 1999). With the 



given Rq and peculiar velocity, the Sun is moving at a veloc- 
ity V Q = (10,226,7) km/s. 

To simulate the motion of DNS in the Galaxy, we randomly 
distribute a population of newly formed DNS in the Galactic 
disk according to a double exponential distribution function 

n(R,z) = n Q exp (j^j ex P (jj^J • ^ 

Here R and z are the cylindrical galactic coordinates, h R and 
h z are the Galactic scale length and height respectively. The 
distribution is normalized to unity by setting no = \ /(Aithjify. 
Since the progenitors of DNS are distributed in the same way 
as the massive stellar binaries in our galaxy, we choose Hr = 
2.8 kpc and h z = 0.07 kpc (see Joshi 2007). For all of the 
observed DNS except PSR B1534+12, we use a population 
of ^10 9 simulated binaries. We use a population of ~10 12 
in the analysis of PSR B1534+12, due to the limiting chance 
of kicking it to the currently high Galactic altitude and si- 
multaneously satisfying the constraints derived from the or- 
bital dynamics of supernova explosions. If the scale height h z 
were larger than 70 pc, it would reduce the difficulty in sat- 
isfying the galactic dynamics constraints of PSR B1534+12. 
A smaller population of simulated binaries would then be 
needed for the analysis of this system. 

The initial center-of-mass velocity of the DNS is obtained 
by summing up the kick velocity Vk, sys imparted to the binary 
center of mass during the second supernova explosion and the 
local Galactic rotational velocity V rot . Given the unknown ori- 
entation of the binary orbit in the Galactic reference frame, 
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the systemic kick velocity V^sys is assumed to be distributed 
isotropically in space. The magnitude of Vk, SV s is drawn from 
the velocity distribution derived in the previous step specific 
to each of the observed systems. 

Starting from the randomly generated initial position and 
velocity, the motion of the DNS is calculated by numerically 
integrating the equations of motion from the DNS birth time 
to the current epoch. The equations of motion are derived 
from the Galactic potential of Carlberg & Innanen (1987) with 
updated model parameters of Kuijken & Gilmore (1989). The 
DNS birth time is chosen according to the DNS age f, gener- 
ated in the first step of the analysis. 

Lastly, we use the measured current position of the ob- 
served DNS to further constraint the ranges of possible Vu and 
Mii. We set the tolerance of this constraint to be within 100 
pc from the observed position. We test the dependence of our 
results on this tolerance level by repeating the analysis with 
a 10 pc, or 50 pc tolerance, but our results do not show any 
significant changes. For PSR J1518+4904 and PSR J1756- 
2251, the position constraint is supplemented with the mea- 
sured current proper motion constraint. For PSR B1534+12, 
PSR B 1 9 1 3+ 1 6, and PSR J0737-3039, the position constraint 
is supplemented with both the measured proper motion and 
NS1 spin tilt angle constraints. 

3. RESULTS 

The analysis outlined in the previous section allows us to 
derive PDFs for the kick velocity Vk, NS2's immediate pre- 
supernova progenitor mass M2;, the spin tilt 8, of NS1, the 
systemic radial velocity V ra d, and the current transverse ve- 
locity Vtrans- Our primary interest is in Vk and Mn, because 
they can shed light on the type of supernova that formed NS2, 
which will be discussed in more detail in the next section. The 
remaining quantities give us some predictions on currently 
unavailable parameters of the 8 known DNS in our galaxy. 
Instead of focusing on the most likely values as studies have 
done in the past, we derive confidence levels for the unknown 
parameters and they are summarized in Table 2. 

Since the Vk and M2; constraints are correlated, we present 
our results in the form of confidence level plots of a 2 dimen- 
sional joint probability distribution of Vj - M 2 ,-, for each DNS 
binary, in addition to the corresponding 1 dimensional PDF's 
of Vk and M He . As pointed out by Willems et al (2006), phys- 
ical parameter constraints are affected by the dimensionality 
of the PDF used, because of inherent correlations between 
parameters and projection effects. This variation is stronger 
when one solely examines the most likely values (i.e., values 
at the peak of the PDF) and weakens as broader confidence 
levels are considered. In what follows, the confidence levels 
of Vk and Mn are derived from the 2D joint PDF. To calculate 
the 2D confidence levels, we first bin the data into a 2D grid 
and normalize the total probability within the 2D grid to unity. 
Then, we add the probability of each bin from the highest to 
the lowest, until the sum best matches the desired confidence 
level. On the other hand, the ranges of V ra d, V trans , and 8 t are 
derived from the corresponding ID PDF. In this case, the con- 
fidence level is found by having the shortest range of bins that 
has a total probability that best matches the desired confidence 
level. 

3.1. PSR Bl 534+1 2 

In this binary, the masses of NS1 and NS2 are 1.33 M and 
1.35 M Q respectively. The orbital period is P h = 0.421 days, 



with an eccentricity e = 0.274. The angle between NSl's spin 
axis and the orbital angular momentum axis is 25.0° ± 3.8° 
or 155.0° ± 3.8°. The binary is currently located out of the 
galactic plane, at I = 19.8° and b = 48.3°, and 1.02 kpc away 
from us. The observed proper motion = 1.3 mas yr" 1 and 
[i-dec. = 25.2 mas yr" 1 . At the measured distance this implies 
V RA . = 6.47 ± 0.32 km s" 1 and V dec . = -121 ± 6 km s _1 , in 
the reference frame of the Sun. Also, we use the spin down 
age, which is 210 Myr given by Arzoumanian et al. 1999, 
as the upper limit of PSR B 1534+12 's age, instead of the 
characteristic age. 

Furthermore, the current orientation of the orbital angular 
momentum axis is also constrained for this sytem. The ori- 
entation axis can be described by the inclination angle / and 
the angle of the orbital ascending node on the plane of the sky 
(fi). The sine of the inclination angle is sin / = 0.975 (Stairs et 
al. 2004). The angle ft is 70° ±20° or 290° ±20°, reck- 
oned north through east. The two solutions correspond to 
cos / < and cos; > respectively (Bogdanov et al. 2002). 
Since tidal effects between the two neutron stars are negligi- 
ble, the orbital angular momentum axis keeps a fixed orien- 
tation in space. From the kinematic analysis of the Galactic 
motion, the constraints on the present-day orbital inclination 
and proper motion can therefore be used to determine the sys- 
temic kick component vf parallel to the post-supernova or- 
bital angular momentum axis. Following Kalogera (1996) and 
Wex et al. (2000), V k \ ys can also be expressed analytically as 



V, 



k,sys 



= y/Ki sintf, = 



^V 2 kz + (Vky+lf 



where 



V r = 



/G(M 2 ,+Mi) 

y A pre SN 

Vkj = ~v ' 

Ml 



Kl 



(M 2i +M X Y 



(2) 



(3) 



(4) 



(5) 



Here, Vkj (j = x,y,z) is the x, y, or z component of the natal 
kick velocity in the frame centered on NS2's progenitor (see 
figure 1 in Kalogera for a graphical representation). Note that 

the sign of vj' needs to be consistent with the sign of Vk z . 

Hence, through equation (2), the V% derived from observa- 
tions gives us a constraint on the y and z components of the 
natal kick, as well as 6, . 

The results are shown in Figures 1 and 2. Right before the 
second supernova, there must have been RLO from the pro- 
genitor of NS2 to NS1, as noted previously by Thorsett et al. 
(2005). As seen in the 2D Vk-Mri joint probability distribution 
in figure 1, the kick imparted to NS2 at birth is between 150 
and 270 km/s, while the progenitor mass M2, is between 1.34 
and 3.40 M , both at 95% confidence. As seen in figure 2, 
polar kicks (i.e., 6' = 0° or 180°) are unlikely, and the pos- 
sible kick directions are asymmetric about the pre-supernova 
orbital plane (i.e., 8' = 90°). This asymmetry comes from the 
constraint on the orientation of the post-supernova orbital an- 
gular momentum axis on the sky. Moreover, we arrive at the 
same conclusion as Thorsett et al. (2005) that 8, ~ 155° is 
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TABLE 2 

Possible Birth Properties and Current Kinematics Parameters of the 8 known DNS in our galaxy 



System 



M 2 ,(M Q r 



y trans (km S )' 



he 



(deg) d 



PSRB1534+12 
PSRB1913+16 
PSR J0737-3039 
PSRJ15 18+4904 
PSR J 1756-2251 
PSRJ1811-1736 
PSRJ1829+2456 
PSR J 1906+0746 



170-260(150-270) 
200 - 410 (190 - 450) 
5 - 50 (5 - 120) 
20-80 (5 - 110) 
5 - 80 (5 - 185) 
0-170(0-310) 
5 - 85 (5 - 225) 
5- 170 (5-510) 



2.00-2.90(1.34-3.40) 
1.40-3.30(1.39-5.00) 
1.25 - 1.55 (1.25 - 1.90) 
1.80-3.30(1.49-4.70) 
1.25 - 1.90(1.25 -2.65) 
1.11 -4.00(1.11 -8.00) 
1.40-2.70(1.36-6.10) 
1.25 -2.90(1.25 -4.80) 



-55 -50 (-110- 145) 
-95 -85 (-255 -315) 
8 - 44 (-20 - 76) 

-30 - 8 (-62 - 48) 
-34 - 32 (-110 -98) 

0- 112 (-90 -210) 
-30 -26 (-126 -66) 
12 - 103 (-176 - 297) 



121 ±6 e 
135 ± 25 e 
23.$*? 
25.4 ± 3.6 e 
6.5 - 28.5 (6.0 - 64.5) 
30 - 1 14 (4 - 202) 
12 - 54 (4 - 102) 
89 - 174 ( 21 - 376) 



25.0 ± 3.8 e 
18 or 162 ± 6 e 
0-2.0(0-7.8) 

2- 18(0-32) 
0-3.4(0- 16.4) 
- 10.2 ( - 57.6 ) 
0-8.5(0-30.0) 
0-6.7(0-51.5) 



NOTE. — The numbers without parenthesis in each entry are at 60% confidence, whereas those with parenthesis are at 95% confidence 
a NS2's immediate pre-supernova progenitor mass 
b Systemic radial velocity 
c Current transverse velocity 
d Spin tilt angle of NS1 

c quantities obtained or derived from observations 
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FIG. 1 . — Kick velocity and progenitor mass of NS2, and present-day kine- 
matic properties of PSR B 1534+12. top: The 2D joint Vt-M 2 , probability 
distribution. The red, yellow and blue colors represent 60%, 90% and 95% 
confidence levels, respectively, middle: The Vj (solid line) and Mij (dashed 
line) PDF. The x-axes of the and Mn PDF are printed on the bottom and the 
top of the plot respectively, bottom: The V tmns (solid line) and V m j (dashed 
line) PDF. The x-axis of the V tmns is printed on the bottom, and x-axis of V raf ; 
PDF is on the top of the plot. 
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FIG. 2.— The 9' PDF 's of PSR B1534+12. 



very unlikely. The current radial velocity of this binary is be- 
tween -110 and 145 km/s, at 95% confidence, and it is likely 
to be moving away from us. Within the observationally con- 
strained range, the V tram PDF is roughly flat. 
When comparing with recent analyses by Willems et al. 

(2004) , Thorsett et al. (2005), and Wang et. al. (2006), 
our limits on Vk and Mx, are consistent wtih Thorsett et al. 

(2005) , and are more constrained than Willlems et al. (2004) 
and Wang et al. (2006), because Willems et al. (2004) did not 
have the spin tilt 6, measurement and Wang et al. (2006) did 
not have the proper motion and position constraints in their 
analysis. Our limits on the angle 6', which is the polar an- 
gle between the natal kick direction and the pre-supernova 
orbital angular momentum, are not as constrained as the lim- 
its shown in Thorsett et al. (2005). Through recent further 
investigation, it has been found that this difference is due to a 
code typo related to the galactic latitude of the source, which 
fortunately affected significantly only the 6' results (Dewey & 
Stairs 2010, private communication). 

3.2. PSRB1913+16 

The masses of NS1 and NS2 in this binary are 1.44 M 
and 1.38 M Q respectively. The orbital period is Pf, = 0.323 
days, with an eccentricity e = 0.617. The angle between the 
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FIG. 3. — Kick velocity and progenitor mass of NS2, and present-day kine- 
matic properties of PSR B1913+16. top: The 2D joint V^-Mn probability 
distribution. The red, yellow and blue colors represent 60%, 90% and 95% 
confidence levels, respectively, middle: The Vj (solid line) and Mij (dashed 
line) PDF. The x-axes of the and PDF are printed on the bottom and the 
top of the plot respectively, bottom: The Vtmns (solid line) and V m j (dashed 
line) PDF. The x-axis of the V, mns is printed on the bottom, and x-axis of V raf ; 
PDF is on the top of the plot. 



NS 1 spin axis and the post supernova orbital angular momen- 
tum axis is 18° ± 6° or 162° ±6°. At the current time, this 
binary is about 300 pc away the galactic plane, at I = 50.0° 
and b = 2.1°, and 8.3 kpc away from us. The observed proper 
motion in right ascension and declination is /ir^ = -3.27 mas 
yr -1 and fi^ec. = -1-04 mas yr -1 . At the measured distance, 
this implies V RA . = -129 ± 25 km s" 1 and V dec , = -41 ± 18 km 
s , in the reference frame of the Sun. In addition, instead of 
the characteristic age, we use the spin down age, which is 80 
Myr (Arzoumanian et al. 1999), as the upper limit on PSR 
B1913+16 's age. 

The results are displayed in Figures 3 and 4. As shown 
in the 2D Vk-Mu joint probability distribution in figure 3, the 
kick imparted to NS2 at birth is between 190 and 450 km s -1 , 
while the progenitor mass M^i is between 1.4 and 5.0 M Q , 
both at 95% confidence. As seen from figure 4, both polar 
and planar kick directions are allowed. The current radial ve- 
locity is between -255 and 315 km s -1 at 95% confidence. 
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FIG. 4.— The 0' PDF 's of PSR B1913+16. 



Within the observationally constrained range, the V, ram PDF 
is roughly flat, with a slight preference for transverse veloci- 
ties of ^120 km/s. Furthermore, although observations show 
a degeneracy in 0, at ~ 18° and ~ 162°, our analysis shows 
that 0, is much more likely to be ~ 18°. 

In addition, although our confidence intervals of Vt and M% 
are consistent with Willems et al. (2004) and Wang et. al. 
(2006), our reported intervals are tighter, because Willems et 
al. (2004) and Wang et al. (2006) reported the allowed ranges 
of Vjt and Mn instead of confidence intervals. 

3.3. PSRJ0737-3039 

The masses of the two neutron stars in this binary are 1.34 
and 1.25 M©, where the more massive one is NS1. They are 
orbiting each other with a period of 0. 102 days, and an eccen- 
tricity of 0.0878. At the current time, the binary is close to 
the galactic plane, at / = 245.2° and b = -4.5°. Through direct 
measurement of geometric parallax, Deller et al. (2009) mea- 
sured the distance to PSR J0737-3039 to be 1.15 kpc, and the 
proper motion in R.A. and declination to be -3.82 and 2.13 
mas yr" 1 respectively. At the measured distance, this implies 
Vra. =-20.8 km s" 1 and Vdec. = 11.6 km s" 1 , in the reference 
frame of the Sun. In addition, Lorimer et al. (2007) con- 
strained the possible age range to be either 70-90 or 170—190 
Myr. We adopt these age ranges in our analysis of PSR J0737- 
3039. 

As the radio pulse profile of PSR J0737-3039A (i.e., NS1) 
has not changed significantly through six years of observation, 
the spin tilt angle of J0737-3039A with respect to the post 
supernova orbital angular momentum axis is believed to be 
small. However, there is no well constrained value for this an- 
gle. Ferdman et al. (2008) and Breton (2008) argue that there 
are two possible models to explain the observed pulse profile: 
the single-cone and two-cone emission model. In this context, 
the number of emission cones means the number of emitting 
cones that intersect with our line of sight. Each model gives 
a different spin tilt angle estimate. Ferdman et al. (2008) es- 
timate this angle to be < 15° for the single-cone model and < 
6.1° for the two-cone model, both at 68.3 % confidence. Due 
to this uncertainty, we take the spin tilt angle 9 t to be < 15° 
in our work, which corresponds to 68.3 % confidence value 
in the single-cone model and 95.4 % confidence value in the 
two-cone model in Ferdman et al (2008). 
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FIG. 5. — Kick velocity and progenitor mass of NS2, and present-day kine- 
matic properties of PSR J0737-3039. top: The 2D joint V k -M 2i probability 
distribution. The red, yellow and blue colors represent 60%, 90% and 95% 
confidence levels, respectively, middle: The Vj (solid line) and Mn (dashed 
line) PDF. The x-axes of the and PDF are printed on the bottom and the 
top of the plot respectively, bottom: The Vtmns (solid line) and V m j (dashed 
line) PDF. The x-axis of the V, mns is printed on the bottom, and x-axis of V raf ; 
PDF is on the top of the plot. 
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FIG. 6.— The orbital separation PDF of PSR J0737-3039 immediately be- 
fore the second supernova. 



PSR J0737-3039 
9,(deg) 

u z 

0.01 o 



0.008 



0.006 



0.004 



0.002 



0.000 

50 1 00 

9' (deg) 

FIG . 7 .— The 0' (solid line) and 9, (dash line) PDF 's of PSR J0737-3039. 
The x-axes of the &' and 9 t PDF are printed on the bottom and the top of the 
plot respectively. 




The results are displayed in Figures 5-7. As shown in the 
2D Vt-M2i joint probability distribution, the kick recieved by 
NS2 (i.e., PSR J0737-3039B) at birth is between 5 and 120 km 
s -1 , while the progenitor mass Mu is between 1.25 and 1.90 
M Q , both at 95% confidence. The probability distribution of 
the kick direction is symmetric about the pre-supernova or- 
bital plane (i.e., 9' = 90°), and planar kicks are more favor- 
able than polar kicks. The PDF of the semi-major axis im- 
mediately before the second supernova (A pre sfi) peaks at two 
narrow ranges: 1.35- 1.40 and 1.45- 1.50 Rq. The first peak 
corresponds to an age of 70-90 Myr, and the second peak to 
an age of 170- 190 Myr. At the current time, PSR J0737-3039 
is more likely to be moving away from us, with a radial ve- 
locity between -20 and 76 km/s with 95% confidence. Within 
the observationally constrained range, the V, rans PDF peaks at 
~20 km/s. Furthermore, RLO from NS2's progenitor to NS1 
before the second supernova explosion cannot be avoided. 

To test the dependence of the presented results on the 
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FIG. 8.— The 6, PDF of PSR J0737-3039, obtained from the same con- 
straints as before except for the upper limit of 8, , which is raised to 60° , to 
test the sensitivity of our results to the adopted 9 t upper limit. 
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adopted upper limit for 6,, we repeated the analysis for the 
upper limits of 5° and 60°. The results do not show any 
significant differences compared to those obtained with the 
15° upper limit. This is not surprising, because the 9, PDF 
as shown in Figure 8, which is obtained from the same con- 
straints as before except raising the upper limit of 6, to 60°, 
peaks strongly at 0° and decays exponentially with increasing 
values of 9, . 

Comparatively, our confidence intervals on Vt and Mn are 
consistent with Stairs et al. (2006) and Willems et al. (2006), 
but are much more constrained than Wang et al. (2006), 
mainly because they did not have the proper motion constraint 
in their analyses. 

3.4. PSRJ1518+4904 

The masses of the two neutron stars in this binary are 0.72 
and 2.00 M©, where the less massive one is NS1. However, 
each one has an uncertainty of ^0.5 M Q . In our analysis, we 
choose Mi to be 1 .23 and Mi to be 1 .49, so that M\ is greater 
than the smallest measured neutron star mass, which is about 
1.2 Mq. The characteristic age is 20 Gyr, which is larger than 
a Hubble time. Thus, we set the upper limit on PSR J1518's 
age to be 10 Gyr instead. The two neutron stars are orbiting 
each other with a period of 8.63 days, and an eccentricity of 
0.249. At the current time, this binary is out of the galatic 
plane, at / = 80.8° and b = 54.3°, and at a distance of 0.625 kpc 
from the Sun. The proper motion in R.A. and dec. is /z«.a. = 
-0.67 and no ec = -8.53 mas yr" 1 . At the measured distance, 
this implies V RA . = -1.99 ± 0.31 and V Dec . = -25.3 ± 3.6 km 
s -1 , in the reference frame of the Sun. There is no spin tilt 
measurement available for this system. 

The results are shown in Figures 9 and 10. As seen in the 
2D Vk-Mii joint probability distribution, the kick imparted to 
NS2 at birth is between 5 and 110 km/s, while the progenitor 
mass Mn is between 1.5 and 4.7 M©, both at 95% confidence. 
Both polar and planar kick directions are possible, but the pla- 
nar kicks have a higher probability. The spin tilt angle 6, is 
smaller than 32° at 95% confidence. The radial velocity is 
found to be between -62 and 48 km/s at 95% confidence, and 
the system is currently likely moving towards us. Within the 
observationally constrained ranges, the V tram PDF is roughly 
flat. 

Our limits of V* and Mu are much more constrained than 
those of Wang et al. (2006), primarily because they did not 
have the proper motion and galactic position constraints in- 
cluded in their analysis. 

3.5. PSR J 1756-2251 

The masses of the two neutron stars in this binary are 1.31 
and 1.26 M©, where the more massive one is NS1. The 
characteristic age is 443 Myr. The two neutron stars are in 
an orbit with a period of 0.320 days, and an eccentricity of 
0.181. At the current time, this binary is close to the galac- 
tic plane, at I = 6.5° and b = 0.95°, and at a distance of 2.5 
kpc away from us. Only proper motion fig a. is measured, 
which is -0.7 mas/yr. At the measured distance, this implies 
Vra. = -8.30 ± 2.37 km/s, in the reference frame of the Sun. 
There is currently no spin tilt measurement available for this 
system. 

The results are shown in Figures 11 and 12. From the 2D 
Vk-Mu joint probability distribution, the kick received by NS2 
at birth is between 5 and 185 km s" 1 , while the progenitor 
mass Mii is between 1.25 and 2.65 M©, both with 95% con- 
fidence. Both polar and planar kicks are allowed, but planar 
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FIG. 9. — Kick velocity and progenitor mass of NS2, and present-day kine- 
matic properties of PSR J1518+4904. top: The 2D joint V k -M 2i probability 
distribution. The red, yellow and blue colors represent 60%, 90% and 95% 
confidence levels, respectively, middle: The Vj (solid line) and Mij (dashed 
line) PDF. The x-axes of the Vj and Mi,- PDF are printed on the bottom and the 
top of the plot respectively, bottom: The Vtmns (solid line) and V m j (dashed 
line) PDF. The x-axis of the V, ra „ s is printed on the bottom, and x-axis of V ra j 
PDF is on the top of the plot. 



kicks are more favorable than polar ones. The spin tilt angle 
of NS1 is smaller than 16.4° with 95% confidence. At the 
current time and at the measured distance, PSR J 1756-2251 
has a total transverse velocity between 6.0 and 64.5 km s -1 , 
and a radial velocity between -110 and 98 km s , both with 
95% confidence. Furthermore, this binary is equally likely to 
be moving towards us as away from us. 

When comparing with Wang et al. (2006), we have tighter 
constrained limits on Vk and Mu, due to the fact that we have 
the additional constraint on the proper motion [Ira.- 

3.6. PSRJ1811-1736 

The two neutron stars have masses of 1.62 and 1.11 M©, 
where the more massive one is NS1. However, each mass 
measurement has an uncertainty of ~0.5 M©. PSR J181 1 has 
a characteristic age of 1830 Myr. The neutron stars are in an 
orbit with a period of 18.8 days, and an eccentricity of 0.828. 
At the current time, this binary is close to the galactic plane, 
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The x-axes of the 0' and 6, PDF are printed on the bottom and the top of the 
plot respectively. 



at I = 12.8° and b = 0.44°, and at a distance of 6.0 kpc away 
from us. Like PSR J 1756-2251, there are currently no spin tilt 
nor proper motion measurements available. 

The results are displayed in Figures 13 and 14. The 2D V^- 
Mn joint probability distribution shows that the kick given to 
NS2 at birth is less than 3 10 km s -1 , while the progenitor mass 
M%i is between 1.11 and 8.0 M©, both with 95% confidence. 
Both polar and planar kicks are allowed, but the planar kicks 
are more favorable than the polar ones. The spin tilt angle 9, 
is smaller than 57.6° with 95% confidence. At the measured 
distance, the total transverse velocity of PSR J181 1-1736 at 
the current epoch is between 4 and 202 km s" 1 , and the radial 
velocity between -90 and 210 km s -1 with 95% confidence, 
which means this binary is likely moving away from us. 

Although our upper boundary of the 95% confidence inter- 
vals of M2i is the same as the upper limit found by Wang et 
al. (2006), the value of our lower boundary is less than their 
lower limit. This is because Wang et al. (2006) put a con- 
servative lower limit of 2.1 M Q on the progenitor mass of a 
neutron star, which we did not. On the other hand, our limits 
on Vk are more constrained than those of Wang et al. (2006), 
as we compute confidence levels instead of allowed ranges of 
solution as Wang et al. (2006). 

3.7. PSRJ1829+2456 

The two neutron stars of this binary have masses of 1 . 14 and 
1.36 M Q , where the less massive one is NS1. However, there 
is an uncertainty of ^0.5 M© in the mass measurements. The 
characteristic age is 12.4 Gyr. For the same reason as PSR 
J 1 5 1 8+4904, we set an the upper limit of 1 Gyr on the age of 
this binary in our analysis. The neutron stars are orbiting each 
other with a period of 1.18 days, and an orbital eccentricity of 
0.139. At current time, this binary is out of the galactic plane, 
at I = 53.3° and b = 15.6°, and it is 1.2 kpc away from us. 
Although there is no spin tilt nor proper motion measurements 
available at current time, Lorimer et al. (2005) derived that the 
transverse velocity is smaller than 118 km s . 

The results are shown in Figure 15. The Vk-Mu joint prob- 
ability distribution shows that the kick given to NS2 at birth 
is between 5 and 225 km s" 1 , while the progenitor mass Mn 
is between 1.4 and 6.1 M© with 95% confidence. The PDFs 
of 0, and 9' are similar to those of PSR J 1756-2251. Both 
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FIG. 1 1 . — Kick velocity and progenitor mass of NS2, and present-day kine- 
matic properties of PSR J1756-2251. top: The 2D joint Vt-Afe probability 
distribution. The red, yellow and blue colors represent 60%, 90% and 95% 
confidence levels, respectively, middle: The Vj (solid line) and Mij (dashed 
line) PDF. The x-axes of the Vj and M?,- PDF are printed on the bottom and the 
top of the plot respectively, bottom: The Vtmns (solid line) and V m j (dashed 
line) PDF. The x-axis of the Vtrans is printed on the bottom, and x-axis of V rac i 
PDF is on the top of the plot. 



polar and planar kicks are allowed, but planar kicks are more 
probable. The spin tilt angle 9, is less than 30.0° with 95% 
confidence. At the known distance, the total transverse veloc- 
ity of PSR J 1829+2456 at the current time is between 4 and 
102 km s -1 , and the radial velocity between -126 and 66 km 
s" 1 with 95% confidence. In addition, the probability that this 
system is moving towards the Sun is ^20% higher than that it 
is moving away from the Sun. 

Comparing our Vk and Mx, results with Wang et al. (2006), 
we have more constrained limits, as they do not include any 
kinematic constraints. 

3.8. PSRJ1906+0746 

The masses of the neutron stars are 1.37 and 1.25 M©, 
where the more massive one is NS 1 . Unlike the other 7 bi- 
naries, the observed pulsar is NS2. This implies the binary is 
relatively young, with a charateristic age of 0.112 Myr only. 
The neutron stars are in an orbit with a period of 0.166 days, 
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FIG. 12.— The 9' (solid line) and 6, (dash line) PDF 's of PSR J1756-2251. 
The x-axes of the 9' and 9, PDF are printed on the bottom and the top of the 
plot respectively. 



and an orbital eccentrirctiy of 0.0853. At current time, this bi- 
nary is close to the galactic plane, at / = 41.6° and b = 0.15°, 
and it is 5.4 kpc away from us. Again, there are no spin tilt 
nor proper motion measurements available yet. 

The results are displayed on Figure 16. The Vk-Mu joint 
probability distribution shows that the kick given to NS2 at 
birth is between 5 and 510 km s" 1 , while the progenitor mass 
M 2i is between 1.25 and 4.80 M Q with 95% confidence. The 
6' and 6, PDF 's are similar to those of PSR J1756-2251. Both 
polar and planar kicks are possible, with the planar kicks hav- 
ing a higher probability. The polar angle 9, is less than 51.5° 
with 95% confidence. At the measured distance, the total 
transverse velocity of PSR J 1906+0746 at the current time is 
between 21 and 376 km s , and the radial velocity is between 
-176 and 297 km s" 1 with 95% confidence. It is most likely 
that this binary is moving away from us. Furthermore, right 
before its supernova, the progenitor of NS2 must have been 
overflowing its Roche lobe, transferring mass to NS 1 . 

Our limits on Vk and M21 are tighter constrained than Wang 
et al. (2006), because they did not have a detailed analysis on 
this system. Instead, they calculated the allowed Vk with an 
assumption that M 2 \ was between 2.1 and 8.0 M Q . 

4. DISCUSSION 

The anisotropic mass and neutrino ejection in a supernova 
explosion gives a recoil natal kick to the newly formed neu- 
tron star. Iron core collapse supernova and ECS produce natal 
kicks of different strength because of the different physical 
conditions of the nascent neutron star immediately before the 
supernovae. 

There are several reviews (Bethe 1990, Mezzacappa 2005, 
Woosley & Bloom 2006, Kotake et al 2006, and Janka et al. 
2007) that summarize our current understanding of the iron 
core collapse supernova. After iron is synthesized in the core, 
the massive star reaches the final stage of hydrostatic nuclear 
burning, because the synthesis of any heavier elements re- 
quires energy instead of releasing energy. When the iron core 
grows by silicon shell burning around the core to a mass above 
the Chandrasekhar mass limit, electron degeneracy pressure 
cannot support the core any more and it starts to collapse. The 
iron core continues to collapse until it reaches nuclear density 
(i.e., pa 10 14 g cm" 3 ). Since nuclear matter has a much lower 
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FIG . 13. — Kick velocity and progenitor mass of NS2, and present-day kine- 
matic properties of PSR J 1 8 1 1-1736. top: The 2D joint Vk-M2i probability 
distribution. The red, yellow and blue colors represent 60%, 90% and 95% 
confidence levels, respectively, middle: The Vj (solid line) and Mij (dashed 
line) PDF. The x-axes of the Vj and Mi,- PDF are printed on the bottom and the 
top of the plot respectively, bottom: The Vtmns (solid line) and V m j (dashed 
line) PDF. The x-axis of the Vtrans is printed on the bottom, and x-axis of V rac i 
PDF is on the top of the plot. 



compressibility, the collapse decelerates and the core bounces 
back because of the increased nuclear matter pressure. This 
drives a shock wave and eventually explodes the outer layers 
of the massive star away in a supernova explosion. 

On the other hand, the ECS scenario is described in the 
reviews of Miyaji et al. (1980) and Nomoto (1987). After 
carbon burning within the core of a post main sequence mas- 
sive star, a O-Ne-Mg core is formed. Neon cannot be ig- 
nited if the mass of the O-Ne-Mg core is less than the crit- 
ical mass of 1.37 M Q for neon ignition. As the O-Ne—Mg 
core is cooled by neutrino emissions, it becomes strongly de- 
generate. When the mass of the core grows by carbon shell 
burning and approaches the Chandrasekhar mass limit, elec- 
tron captures onto 24 Mg and 24 Na take place. This leads to 
a decrease in the electron degeneracy pressure and the Chan- 
drasekhar mass limit, which induces a rapid contraction of the 
core. The rapid core contraction ignites the oxygen deflagra- 
tion, which incinerates materials into nuclear statistical equi- 
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librium (NSE). Due to the rapid electron capture onto NSE 
elements, the collapse of the core accelerates. Similar to the 
iron core collapse picture, when the core reaches the nuclear 
matter density, it bounces back and drives a shock wave that 
eventually explodes the star in a supernova. Gutierrez et al. 
(2005) show that the abundance of 24 Mg in O-Ne-Mg core 
needs to be greater than 15% , in order to have NSE devel- 
oped through electron captures. However, in simulations of 
evolving massive AGB stars by Siess (2007), the abundance 
of 24 Mg after carbon burning is smaller by a factor of ~10 
than what is required to drive an explosion by electron cap- 
tures. Therefore, the mechanism of developing an ECS is not 
a settled issue yet. In addition, Podsiadlowski et al. (2004) 
studied the binary evolution and the dynamics of core col- 
lapse, and suggested that ECS could only occur in interacting 
binaries, but not in single stars. 

In recent simulations of supernova explosions from the col- 
lapse of O-Ne-Mg cores by Kitaura et al. (2006), they 
found these supernovae are powered by neutrino heating and 
neutrino-driven wind of the nascent neutron star. Scheck et 
al. (2006) showed that the shock wave in these supernovae 
can propagate outwards on a relatively short timescale after 
the rebound of the core, which means the non-radial hydro- 
dynamics instabilities do not have time to merge and grow to 
global asymmetry before the anisotropic pattern freezes out in 
the accelerating outward motion of the shock wave. As a re- 
sult, the natal kick of an ECS due to anisotropic mass ejection 
is expected to be fairly small (see Scheck et al. 2004). Fur- 
thermore, as the O-Ne-Mg core that eventually explodes in 
an ECS is not massive enough to ignite neon, the mass of the 
neutron star progenitor should be less than the progenitors 
that explode in an iron core collapse supernova. 

The present analysis is the first to include all currently 
known DNS systems and account of all orbital and kinematic 
constraints. We also employ a novel method for dealing with 
the uncertainty due to the un-measurable radial velocities and 
we focus on the derived PDFs and associated confidence lev- 
els, instead of just the most likely values, which can be mis- 
leading at times. Our results are summarized in Table 2; the 
derived constraints are consistent with earlier studies (when 
available) but typically limit parameters to narrower ranges, 
and in this sense they represent the best available constraints 
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FIG. 15. — Kick velocity and progenitor mass of NS2, and present-day kine- 
matic properties of PSR J1829+2456. top: The 2D joint V k -M 2i probability 
distribution. The red, yellow and blue colors represent 60%, 90% and 95% 
confidence levels, respectively, middle: The Vj (solid line) and Moi (dashed 
line) PDF. The x-axes of the Vj and Mi, PDF are printed on the bottom and the 
top of the plot respectively, bottom: The Vtmns (solid line) and V m j (dashed 
line) PDF. The x-axis of the Vtrans is printed on the bottom, and x-axis of V ra j 
PDF is on the top of the plot. 



on the formation of neutron stars in DNS systems. 

In the context of our current understanding of massive star 
core collapse, we can use our results to draw a number of con- 
clusions: 

(1) PSR J0737-3039 has a V k upper limit of 120 km/s and a 
Mu upper limit of 1.9 M Q at 95% confidence. Therefore, the 
formation of NS2 (i.e., pulsar B) likely occurred through an 
ECS event. This is consistent with the speculation of Podsi- 
adlowski et al. (2005) based on the space velocity and orbital 
eccentricity of PSR J0737-3039. For PSR B 1534+12 and PSR 
B1913+16, V k are 150 - 270 km/s and 190 - 450 km/s (95%) 
respectively, which means NS2 in both system must have re- 
ceived a significant recoil natal kick at birth. PSR B1534+12 
has a M 2i upper limit of 3.4 M Q (95%), while PSR B 19 13+ 16 
has a Mri upper limit of 5.0 M2, (95%). Because of the rela- 
tively high Vk ranges and Mn upper limits, NS2 in both sys- 
tems are probably formed through an iron core collapse su- 
pernova event. 
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FIG. 16. — Kick velocity and progenitor mass of NS2, and present-day kine- 
matic properties of PSR J 1906+0746. top: The 2D joint V k -M 2 i probability 
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confidence levels, respectively, middle: The Vj (solid line) and Mij (dashed 
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(2) PSR J 15 18+4904 has a V K upper limit of 1 10 km/s and a 
Mn upper limit of 4.7 M Q at 95% confidence. Even though 
NS2 likely received a low recoil velocity at birth, we cannot 
firmly conclude which type of supernova occurred during the 
formation of NS2, due to the relatively high Mji upper limit. 

(3) PSR J1756-225 1 has a V k upper limit of 80 km/s and a Ma- 
upper limit of 1 .9OM at 60% confidence, hence the forma- 
tion of NS2 might possibly relate to an ECS event. However, 
the Vk and Mr, upper limits are 185 km/s and 2.65 M Q re- 
spectively at 95% confidence, which means the likelihood of 
NS2 formed through an ECS event is lower than that of PSR 
J0737-3039. 

(4) PSR J 1 8 1 1 - 1 736, PSR J 1 829+2456, and PSR J 1 906+0746 
have Vk upper limits of several hundreds of km/s, and Mn 
upper limits ranging from 4.80 to 8.0 M Q (i.e., the conven- 
tional limit on neutron star progenitor mass) at 95% confi- 
dence. Since both low and high V* and My are possible, we 
cannot conclude which type of supernova event is more favor- 
able for the formation history of NS2 in these systems. 

(5) For all of the known DNS except PSR B 1534+12, all avail- 
able constraints are consistent with imparting a polar or planar 
kick to NS2 at birth. For PSR B1534+12, the kick direction is 
constrained to be polar. 

(6) For PSR B1534+12, PSR J0737-3039 and PSR 
J 1906+0746, the pre-supernova orbit is so tight that they can- 
not avoid a RLO from NS2 progenitor to NS1 just before the 
second supernova explosion. 

(7) Furthermore, despite the low orbital eccentricity in PSR 
J 1829+2456 and PSR J 1906+746, NS2 could have a high 
progenitor mass and have received a high recoil natal kick 
at birth, as shown in the 2D joint PDF of these systems (see 
Figures 15, and 16). 

(8) We also tested the dependence of our PSR J0737-3039 
results on the current upper limit of the spin-orbital misalign- 
ment angle 8, of NS1 (i.e., pulsar A). As shown in Figure 8, 
the 9, PDF peaks strongly at ~ 0°, so choosing an upper limit 
of 15° or 60° does not affect our results noticably. 
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